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ABSTRACT 

We present the results of wide-field deep JHK imaging of the SSA22 field using MOIRCS instrument 
equipped with Subaru telescope. The observed field is 112 arcmin 2 in area, which covers the z = 3.1 
protocluster characterized by the overdensities of Lya emitters (LAEs) and Lya Blobs (LABs). The 
5 a limiting magnitude is Kab — 24.3. We extract the potential protocluster members from the 
A-selected sample by using the multi-band photometric-redshift selection as well as the simple color 
cut for distant red galaxies (DRGs; J — Kab > 1-4). The surface number density of DRGs in our 
observed fields shows clear excess compared with those in the blank fields, and the location of the 
densest area whose projected overdensity is twice the average coincides with the large-scale density 
peak of LAEs. We also found that AT-band counterparts with z p hot — 3.1 are detected for 75 % 
(15/20) of the LABs within their Lya halo, and the 40 % (8/20) of LABs have multiple components, 
which gives a direct evidence of the hierarchical multiple merging in galaxy formation. The stellar 
mass of LABs correlates with their luminosity, isophotal area, and the Lya velocity widths, implying 
that the physical scale and the dynamical motion of Lya emission are closely related to their previous 
star-formation activities. Highly dust-obscured galaxies such as hyper extremely red objects (HEROs; 
J — Kab > 2.1) and plausible AT-band counterparts of submillimeter sources are also populated in 
the high density region. 

Subject headings: galaxies: formation — galaxies: high-redshift — galaxies: evolution — cosmol- 
ogy:observations — galaxiesxlusters: general 



1. INTRODUCTION 

Since the late 1990's, overdense regions of high redshift 
galaxies at z > 2 have been discovered and extensively in- 
vestigated by many authors. The most notable examples 
are high density regions of Lya emitters (LAEs), Lyman 
break galaxies (LBGs) and their resemblances over the 
redshift range of 2 - 5 (e.g. Steidel et al. 1998, 2000; 
Shimasaku et al. 2003; Ouchi et al. 2005). These stud- 
ies revealed the properties of the on-going star formation 
activities in the high density regions of the universe at 
high redshift. However, the galaxies selected by the rest- 
frame ultra-violet (UV) emission may not be the only 
representative population in such dense regions. For ex- 
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ample, sum of the stellar mass of a hundred of LAEs is 
only 10 10-11 M©, as the typical LAEs are known to be 
less massive objects with ~ 10 8-9 M (Gawiser et al. 
2007; Ono et al. 2010), which corresponds to only one 
typical massive galaxy. 

On the other hand, the overdensities around high-z ra- 
dio galaxies and QSOs, which are expected to evolve into 
giant elliptical galaxies, are also observed so far (Kurk et 
al. 2004; Venemans et al. 2007; Overzier et al. 2009; 
Hatch et al. 2011). From recent near- to mid- infrared 
and submillimeter observations, the high-density regions 
of relatively mature galaxies or highly-obscured dusty 
starbursts are indeed detected around the radio galax- 
ies, or X-ray selected cluster candidates (Kodama et al. 
2007; Daddi et al. 2009; Gobat et al. 2011; Tanaka et 
al. 2011). The samples are still very limited, however, 
and the general properties of the protoclusters at high 
redshift are far from understood. 

The protocluster at z — 3.09 in and around the SSA22 
field was first discovered by Steidel et al. (1998, 2000) as 
a significant redshift peak of LBGs with the overdensity 
of <5gai ~ 5.0 over the 20 Mpc-scale comoving volume. 
Later narrow-band observations have found the large- 
scale filamentary structure of LAEs extending over 60 
Mpc (the SSA22-sbl field; Hayashino et al. 2004). Ya- 
mada et al. (2012) revealed that the overdensity of LAEs 
is in fact the most prominent structure in the larger sur- 
vey area of 1.38 deg 2 , corresponding to the comoving 
volume of 10 6 Mpc 3 . 

In this field, along the structure of LAEs, there also 
locate the extremely large, giant Lya clouds, Lya-blobl 
(hereafter LABI) and blob2 (LAB2) (Steidel et al. 2000), 
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and the other 33 LABs (Matsuda et al. 2004, 2007). The 
objects are considered to be closely related to the early 
phase of massive galaxy formation in the overdense re- 
gions at high redshift (e.g., Yamada 2009). LABs typ- 
ically have the physical extent of 30 — 150 kpc and the 
Lya luminosity of > 10 erg s _1 . They are mostly found 
in the high-density regions of star-forming galaxies at 
z = 2 - 5 (Matsuda et al. 2004, 2011; Saito et al. 2008; 
Yang et al. 2010), while the most distant LAB is discov- 
ered by Ouchi et al. (2009) at Z = 6.6. The origin of the 
strong and extended Lya emission is still unclear; it may 
be powered by cold gas accretion (Haiman et al. 200; Go- 
erdt et al. 2010a), or by galactic superwinds (Taniguchi 
& Shioya 2000; Mori & Umemura 2006), or by photoion- 
ization induced by intense starbursts or active galactic 
nuclei (AGN) (Chapman et al. 2001; Geach et al. 2009). 
Recent studies based on X-ray and mid-infrared observa- 
tions suggests that a significant amount of LABs has an 
indication of AGN which can power the extended Lya 
emission (Lehmer et al. 2009; Geach et al. 2009; Webb 
e t al. 2009). 

[Uchimoto et al. ( 2008 ) carried out a deep near-infrared 
observati on in the southern part of the SSA22a field 
( Steidel et al.||2000 ). In the field (hereafter referred as 
SSA22-M1), they found density excesses of distant red 
galaxies (DRGs; Franx et al. 2003) and the photometric 
redshift (hereafter photo-z) selected objects in the vicin- 
ity of LABI and LAB2. Moreover, they showed that the 
stellar mass of K-selected objects associated with 8 LABs 
are in the range of M* = 10 9 — 10 11 M , and they are 
correlated with their luminosity of Lya emission. Since 
the observed field of view is not large enough to cover the 
entire high-density region of LAEs, however, the stellar 
mass distribution over the entire high-density structure 
was not clear very much. 

In this paper, we present the results of the new deep 
and wide near-infrared (NIR) imaging of the protocluster 
at z = 3.09 in the SSA22-sbl field by using Multi-Object 
Infra-Red Camera and Spectrograph (MOIRCS; Suzuki et 
al. 2008) mounted on 8.2 m Subaru Telescope. We here 
construct the large NIR-selected sample of protocluster 
galaxies down to K = 24.3 along the filamentary, high- 
density structure of LAEs at z = 3.09. The limiting 
magnitude corresponds to the stellar mass > 3x 10 9 M Q 
assuming the typical SED of LBGs at z ~ 3. We will 
focus on the stellar mass assembly in the protocluster at 
z = 3.1 in this paper. Our scientific goal is to reveal 
how much of the stellar mass has already formed and 
assembled in the protocluster. In other words, we would 
like to see the relation between the structure traced by 
the star-formation activity and that traced by the stellar 
mass. 

We will describe the observation and the data analysis 
method in section [2] The sample selection by the photo- 
metric redshift (photo-z) as well as by the simple color 
cut will be described in section [3] The properties of all 
the selected objects and AVbana counterparts of LABs 
will be explained in sectionHJ In section[5]we will discuss 
our results in terms of the mass assembly history in the 
protocluster. 

We use the cosmological parameter values f_M =0.3, 
S1 A = 0.7, and H o =70 km s _1 Mpc - 1 throughout 
this paper. All the magnitud e values are in AB sys- 
( |Oke fc Gunn~~||1983| [Fukugita et~aLl|1996| ), un- 
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Fig. 1.— The observed area in the SSA22-Sbl field. The rect- 
angles (red) indicate each field of view, in which the names of the 
observed fields are also shown. Contours (black) are la, 1.5(7, 2(7 
density levels of LAEs. LABs are indicated with circles (blue). The 
cross shows the density peak of LAEs identified in Yamada et al. 
(2012). 

less explicitly noted. Conversions between the Vega and 
the AB systems are made using the following equations: 

J = Jvoga + 0.95, H = ffvcga + 1-39, K = Avcga + 1-85. 

2. OBSERVATIONS AND DATA 

J, H, and A s -band images were obtained by using 
MOIRCS equipped with Subaru telescope, whose field of 
view is 4' x7' . The observe d fields are located in SSA22- 
Sbl (Matsuda et al. 2004). The observations were car- 
ried out on June and August 2005, July 2006, September 
2007, and May 2008. The summary of the observations 
is listed in Table [T] 

We observed the six regions, which are referred as 
SSA22-M1 to M6. The schematic feature is shown in 
FigJT} The name of the observed field is also shown in 
the figure. The area of the three fields, M3, M4, and M5, 
are 4' x 3.5', since the half of the FOV was not available 
in 2007 due to the detector trouble. The early results 
for the SSA22-M 1 field, correspondi ng to the southern 
part of SSA22a ( Steidel et al. 1998), was published in 
Uchimoto et al. ( 2008[ ). The additional data taken for 
the SSA22-M1 fields also analyzed in this paper. 

The images are reduced in standard manner wi th the 



MCSRED so ftware package PI (Tanaka I. et al.). Uchi- 



|moto et al.| |2008) can be referred for the detailed de 
scription. The effective area of the final image is 111.8 
arcmin 2 . The typical image size (FWHM) ranges from 
0".4 to 0".6, and the limiting magnitude ranges from 24 
to 25 mag. The details are listed in Table [2j 
For the detection an d photometry, we us e SExtrac- 



1996). The ob- 



tcm 



tor (version 2.3) ( jBertin fe Arnouts 
jects which have more than 16 connected pixels above 
1.5 a in surface brightness are selected. The SExtrac- 
tor MAG_AUTO value is adopted as the A-band pseudo 

11 http://www.naoj.org/staff/ichi/MCSRED/mcsred_e.html 
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total magnitude. We calibr ated our NIR data to the 
UKIRT photometric system ( |Tokunaga et al.||2002[ ). 

For the J — K color measurement, all the images of 
the SSA22-M1, M2, M3, and M6 fields are smoothed to 
be matched with 0".53, the largest PSF of the J-band 
images in Ml. On the other hand, FWHM in M4 and 
M5 are slightly larger as shown in Table [T] We therefore 
matched the PSFs in M4 and M5 to 0".67. The J - K 
color was measured by a l.l"(9-pixels)-diameter aper- 
ture in Ml, M2, M3 and M6, and by a 1.4"(12-pixels)- 
diameter aperture in M4 and M5, and we applied a small 
aperture correction less than 0.1 mag to the aperture 
magnitudes of the sources in M4 and M5. Fig[2| shows 
the J — K vs. K color-magnitude diagrams for the four 
regions with better image quality, and FigJ3] shows that 
in the two regions with relatively poor image quality. In 
SSA22-M1, M2, M3 and M6, the J-band images are deep 
enough to select the DRG sample (J — K > 1.4) down to 
K = 24.0. The completeness in the shallowest M4 field 
is still 90 % at K = 23.8. 

To estimate the photometric redshift, we also use the 



optical BV Ri' z' -bands data obtained by Matsuda et al.| 
(2004) with Subaru/Suprime-Cam as well as w*-band 
dalataken by CFHT/MegaCam (taken by P.I. L.Cowie). 
As the image size in optical and ultraviolet (UV) data is 
1".0 in FWHM, we smoothed the J, H, and K band im- 
ages to be matched with the optical ones and measured 
the colors in 2" aperture for the purpose. 

The Galactic extinction was corrected by adopting the 
av erage value at the field, E(B — V) — 0.08 (Schlcgel ct 
aI] [T998l ). 



3. SAMPLE SELECTION 

We here select the potential protocluster members 
from the /^-selected sample. Since if-band corresponds 
to rest-frame U-band at z ~ 3, the if -band lu minosity of 
z = 3 galax ies nearly reflect their stellar mass ( Kajisawa 
et al.|2009"j) . The nominal stellar mass limit at K = 24.0 



corresponds to ~ 10 95 M Q . 

We adopted the two kinds of selection criteria; 
one is the photometric redshift analysis by using 
u* BV Ri 1 z' J H K data, and the other is the simple color- 
cut method. While full photometric redshift analysis is 
useful, it is still inevitable to have some scatters and 
a fraction of catastrophic failure if the photometric er- 
rors are ineligible and the limited model templates are 
used. As the photometric redshift needs the optical im- 
ages with poorer image quality, it suffers from the rela- 
tively large photometric error due the larger photometric 
apertures (see above). A simple color selection such as 
DRGs thus works as a complementary method not only 
as they are simple and reproducible but also as the pho- 
tometric error is relatively small especially for the fainter 
objects. 

3.1. Photometric Redshift of the K-selected Sources 

Using the multi-band photometric data, we obtained 
the photome tric redshift of the if-s elected sources using 
hyperz code ( |Bolzonella et al.||2000[ ). The SED fitting is 
performed with the redshift, spectral type, age, and dust 
extinction as the free parameters. The best fit SED is 
determined by minimizing \ 2 value. The template spec- 
tra we used are those for the simple stellar population 



with the fixed (solar) me tallicitiy from GALAXEV by 
Bruzual fc Charlot| p03] ) . 

Fig[3] shows the photometric redshifts for the spectro- 
scopic redshifts obtained by the previous literature and 
the NASA/IPAC Extragalactic Database (NED). The 
relative errors, (z p hoto — z sp)/(l + z sp ), are also shown 
in the bottom panel. The errors arc similar if compared 
with those in the previous studies (e.g. Ichikawa et al. 
2007). The uncertainty of Az ~ 0.5 at z ~ 3 remains 
since the Balmer/4000 A break is shifted in the middle 
of ii-band. 

The median photometric redshift of 26 LBGs with 
3.06 < z S pcc < 3.13 is Zphot = 3.11 and the standard 
deviation is 0.08. 92 % of the LBGs is in the range of 
2.6 < Zphot < 3.6. We therefore pick up the objects with 
2.6 < Zphot < 3.6 as the candidate galaxies at z = 3.1. 
We finally selected 474 objects brighter than if =24. 



3.2. Selection of Distant Red Galaxies (DRGs) 

The red J — K color of DRGs is due to the Balmcr 
or 4000A break of g alaxies at z > 2, or to the dust ex- 



tinction dFranx et al. 2003; Reddy et al. 2005 Forster 
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|20TJ4| . JWuyts et al |J 2~009 confirmed by 



spectroscopic observation that DRGs with Ky cga < 22.5 
(-Kab < 24.35) are dominated by z > 2 galaxies. In the 
GOODS-North field, 86 % of DRGs are in the range of 
2 < Zphot < 4 (Kajisawa et al. 2011), and the median 
photometric redshift is z p hot = 2.5. 



We used the same criteria as in Kajisawa et al. (2006), 
which corresponds to J — Kab > 1.4 in our photometric 
system. In our data, we sampled 356 DRGs with K < 
24 in the whole observed area. The distribution of the 
photometric redshifts for DRGs is shown in Fig[5] We 
find that 70 % of DRGs in the SSA22 area are located 
at 2 < Zphot < 4. The median photometric redshift is 
Zphot =2.6. 

We compared the current sample of DRGs with the 
previously known sample of LBGs (Steidel et al. 2003). 
In the area of SSA22a, which has a substantial overlap 
with our MOIRCS fields, they listed in total of 171 LBGs 
including those without spectroscopic redshift (29 ob- 
jects have the redshift between z=3.06 and 3.12). Of 
them, only the four objects (SSA22a-C54, SSA22a-M14, 
SSA22-MD39, and SSA22a-aug96C20) meets the J - K 
color criteria. The if-band counterpart of SSA22a-M14a 
is, however, displaced from the optical source and may 
not be the same object or the same region in a galaxy. 
The spectroscopic redshift of SSA22a-aug96C20 is 1.357. 
No spectroscopic redshift is available for SSA22a-C54 
and SSA22a-MD39. 

We also checked the optical colors of our sample of 
DRGs to see how many have the colors match the LBG 
criteria. Among the 356 objects with K < 24, 7 ob- 
jects with R < 25.5, magnitude limit for the spectro- 
scopic sample of Steidel et al. (2003), and 23 objects 
with R > 25.5 satisfies the LBG criteria adopted here, 
u* - V > 2.9(7 - R) + 0.37, and V - R < 0.5, and 
u* — V > 1.4. These numbers should be referred care- 
fully as the color selection only properly works for those 
objects with high photometric accuracy (> 10cr in V & 
R). The optical colors of those optically faint red objects 
are largely scattered due to the photometric errors. In 
conclusion, we found that the overlap between the cur- 
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Fig. 2. — J — K vs. K color-magnitude diagram of the if-selected objects in SSA-M1, M2, M3, and M6. The diameter of aperture 
photometry is 9 pixels (l".l) is used for color measurement. The objects fainter than the J-band detection limit is show by the arrows. 



rent samples of DRG and LBG is very small. 

3.3. Selection of Hyper Extremely Red Objects 
(HEROs) 

We also sampled hyper extremely red objects (HEROs; 
Totani et al. 2001) with J — Kab > 2.1 as a subset of 
DRGs. Totani et al.(2001) suggested that the extremely 
red color in J — K is well explained by primordial ellip- 
tical galaxies reddened by dust and still in the starburst 
phase at z ~ 3. Indeed, in the GOODS-North field, all 
of the HEROs are in the range of 2 < z p hot < 4 (Kaji- 
sawa et al. 2011), and the median photometric redshift 
of HEROs is Zp hot = 3.0. 

The number of HEROs in the observed field is 31. We 
find that 94% (29/31) of HEROs in our field are classified 



as 2 < Zphot < 4. The median photometric redshift of 
HEROs is Zp ho t = 2.7. 

It should be noted that the result of the SED fitting for 
the HEROs tends to include relatively large uncertainty 
due to their extreme red color. 55% of the HEROs are 
not detected in J-band at 2 a level (see Figs J2] and [3|. 
The number of the HEROs which are detectea only in 
K-band is four while the rest of the objects are detected 
at least in the relatively deeper i?-band image. 

4. RESULTS 

4.1. The Distributions of the K-selected Objects in the 
SSA22 field 

The distributions of the if-selected objects with 
Zphot = 2.6—3.6 and DRGs on the sky are shown in FigM 
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The same as Fig[2] but for the X-selected objects in SSA22-M4 and M5. The diameter of aperture photometry is 12 pixels 
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Fig. 4. — Comparison between photometric redshift and spec- 
troscopic redshift for the if-selected objects in the SSA22-M1, M2, 
M3, M4, M5, and M6 fields. The relative errors (z pho to —z ap )/(l + 
2 S p) are indicated in the bottom panel. 



Fig. 5. — Histogram of the estimated photometric redshift of 
DRGs and HEROs with K < 24.0. The red histogram shows the 
numbers of DRGs and the hashed histogram shows the numbers of 
HEROs. 



where LABs (Matsuda et al. 2004) and LAEs (Hayashino 
et al. 2004) are also plotted. The figure shows the sam- 
ples with the magnitude limit of K < 24.0. The dotted 
lines are the region observed with MOIRCS. We note 
that the noisy area at the center of the SSA22-M6 field 
is removed for the DRGs selection to avoid the spurious 
detection. The solid contours delineate the high-density 
region of the LAEs in Hayashino et al. (2004). 

In order to investigate whether there is an excess in 
number density of the galaxies in the SSA22 protoclus- 



ter compared with the blank fields, we first investigate 
the surface number density of DRGs. The cumulative 
and differential number counts of DRGs in the whole 
field are shown in FigjTJ The DRG counts show the ex- 
cess at K > 23 compared to those in all the three blank 
fields, namely the GOODS-North field (MOIRCS Deep 
Survey; Kajisawa et al. 2006, 2011), the GOODS- South 



field, and the HDF-South field ( jGrazianet al.|2006[ ). The 
completeness fraction and the 90 % limits in SSA22-M1 
(deepest) and M4 (shallowest) are also plotted. The cu- 
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Fig. 6. — (upper) The distribution of photo-z selected galaxies (z D hot = 2.6 — 3.6) (open pentagons, red), (bottom) The distribution 
of DRGs (open circles, red). LABs (Matsuda et al. 2004) and LAEs (jHayashino et al.|2004| are indicated with large open circles (blue) 
and open squares (cyan), respec tively. The observed fi elds with MOLKUS are surrounded with the dotted lines. The contour shows the 
high-density region of the LAEs JHayashino et al.|2004[l. The cross is the density peak of LAEs identified in Yamada et al. (2012). 
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Fig. 7.— (left) Cumulative number counts of DRGs in SSA22-MOIRCS (Ml, M2, M3, M4, M5, and M6) (open circles) in comparison 
with those for DRGs in GOODS-N (crosses), GOODS-S (open squares), and HDF-S (open triangles). The detection completeness in 
i\"-band in SSA22-M1 (deepest) and M4 (shallowest) is shown with the smaller dotted line. The Poisson errors are shown by the error bars, 
(right) Differential number counts of DRGs in SSA22-MOIRCS. 



initiative number density of DRGs in the entire observed 
field is 3.1 arcmin 2 down to K = 23.85, which is ~ 33 % 
larger than that in GOODS-N. 

FigjH] shows the DRG number counts in each field of 
view. The number density in GOODS-N and the 90% 
completeness fraction are also plotted in the figure. The 
densest fields are SSA22-M3 and M6, which are over- 
lapped with each other, followed by SSA22-M4, M2, M5, 
and Ml. The number density of the DRGs in M6 is 
2.1-1.8 times larger than that in GOODS-N at K=23- 



In addition, the excess of HEROs is found in this 

Our observed fi elds cover 
AE s obse rved by Matsuda 



23.5 

field (see section |5.1| and [5.3 ) 
the highes t densit y peak oF 



et al. (|2011|) and Yamada et al.| (|2012) in the 1.38 deg z 

(Fig tt The outstanding density peak of 

3 arcminutcs lies in SSA22-M3 



survey area 

LAEs with the size of 

and M6, whose location is shown by the cross symbol 

in FigsfTJ anc|6j The fact that the largest overdensity 

of DRGs appears at around the highest density peak of 

LAEs suggests that the stellar mass distribution traces 
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the structure of LAEs in the scale of several arcminutes, 
or a few Mpc at z = 3.1. 

We also investigated the surface number density of 
photo-z selected objects and found that the counts of 
the objects at 2.6 < z p hot < 3.6 in the SSA22 fields are 
larger than those in the GOODS-N field (see Fig. 9 be- 
low). The direct comparison of the two samples should 
be discussed carefully, however, since the depth, or the 
photometric uncertainty of them is different. 

These results for the if -selected galaxies can be com- 
pared with the surface density of LBGs in the SSA22 
fields. Steidel et al. (2003) found the average surface den- 
sity of the color-selected LBG candidates with R < 25.5 
to be 1.72 arcmin -2 . While the surface density of the 
11 fields observed with the Palomar 200-inch telescope 
ranges from 1.10 to 2.15 arcmin -2 , the density in SSA22a 
and SSA22b are 1.88 and 1.15 arcmin -2 , respectively. 
The higher contrast of the protocluster in the if-selected 
galaxies is due to the fact that the more stellar massive 
galaxies traced by the if -band sample have intrinsically 
much stronger clustering properties and may be easier to 
be dominated by a single large structure. 

4.2. K-band Counterparts of Lya Blobs 

As seen in Figj6j photo-z selected objects and DRGs 
are frequently seen around the position of LABs. We 
investigate the spatial correlation between these K- 
selected objects and rest UV-selected galaxies (LABs and 
LAEs) in the field. Fig(9] shows the surface number den- 
sities of if-selected objects as a function of a radius from 
LABs and LAEs. As shown in the figure, we found that 
the if-selected objects are correlated with LABs within a 
radius of 6", which corresponds to 50 kpc scale at z = 3.1. 
This implies that these objects are likely to associate with 
LABs. On the other hand, no significant correlation be- 
tween LABs and if-selected objects at the larger scale is 
detected. 

In contrast to LABs, LAEs are not correlated with K- 
selected objects in scale of less than one arcminute while 
the large-scale density peak s of LAEs and if-selected ob- 
jects coincide ( section J4.1[). Most of the LAEs are not 



detected in the if-selecte d ima ges to the limit except for 
some HEROs (see section 5.3.1 ). The expected typical lu- 



minosity of LAEs in the SSA22 protocluster is K ~ 26.3 
which is estimated from the stacked image of LAEs that 
are not detected in if-band individually. The result is 
consistent with that of LAEs in a field environment (Ono 
et al. 2010). 

We here describe the detailed near-infrared properties 
of LABs. The number of the LABs (Matsuda et al. 2004) 
in our observed fields is 20, and the if-band co unt erpa rts 
are detected for 75 % (15/20) of them. Figs [lO] and |TI] 
show the images of the 15 LABs, which have plausible LC 
band counterparts. Each panel has a length of 23. "4 on a 
side, which corresponds to 180 kpc at z = 3.1. Circles are 
photo-z selected objects and DRGs with z p hot = 2.6 — 3.6 
are colored red. Red squares are DRGs with z p hot < 2.6 
or Zphot > 3.6. The contours are isophotal area of Lya 
emission with surface brightness of 28.0 mag arcmin -2 
obtained by the narrow-band observation in Matsuda et 
al. (2004). 

Since the expected average surface number of DRGs 
or photo-z selected objects is only ~ 0.03 in the median 
Lya isophotal area of 23 arcsec 2 , or 0.0012 arcsec -2 in 
surface density, it is interesting that a significant number 
of if-selected objects are detected around the LABs in 
our field. Table[3|shows if -band magnitude, J—K color, 



^•ph ot 



and 



11 



'spec 



of the if-selected objects in Figs 10 and 



Spectroscopic redshifts of 11 photo-z selected objects 
including 9 LBGs listed in Tabic [3] have been confirmed 
to be z sp cc = 3.1. 

On the other hand, we detect no if -band counterparts 
for 5 LABs, LAB9, 19, 25, 26, and 35, as shown in FigJTJJ 
Please note that we plot all photo-z selected objects and 
DRGs even below the nominal magnitude limit. Detailed 
description for each object is given below. 

4.2.1. LABs with K-band Counterparts 

We identified the if-band counterparts of the LABs, 
which are numbered in FigsjlO] and [TT] There are the 
two types of LABs among them, those with a single or 
multiple if -band counterparts. LAB3, 5, 8, 14, 20, 24, 
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Fig. 9. — (upper) The surface number density of A"-selected 
objects with z p h ot = 2.6 — 3.6 in a circular ring of 0.1 arcmin width 
as a function of radial distance from LABs or LAEs. The surface 
number densities from LABI and LAB2, LABs except LABI and 
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open squares (green) , respectively. The mean densities of DRGs in 
SSA22 and GOODS-N are indicated with dotted and solid lines, 
respectively, (bottom) The surface number density of DRGs in a 
circular ring of 0.1 arcmin width. 



and 31 are associated with a single component. The other 
8 LABs, LABI, 2, 7, 11, 12, 16, 27, and 30, have the 
multiple components. 

LABI has six if-band components within the Lya 
halo. #4, #5 and #6 are newly detected in this deeper 
.ftT-band image. While the best-fit SED of #5 gives the 
photo-z outside the criteria z p hot = 2.6 — 3.6, it is still 
a possible counterpart at z ~ 3.1 if the large photomet- 
ric errors in NIR are taken into account. There is a 
DRG (#2) near the center of the LAB, and one LBG at 
z = 3.109 (#3) in the southwest from the center. As sug- 



gested in U08, a deficit of Lya emission is fo und at the 
position of #1. No X-ray source is detected (Lehmer et 
al. 2009), and there is no indication of AGN at the mo- 
ment. MIPS 24 /im emission is detected at the position 
of #1 as well as the region in the east from the center. 
The spectroscopic redshift at the Lya e mission peak and 
at th e position of #2 agree to be 3.097 (Ohyama et al. 
20031) . Near #4, the optical conti nuum object is detectecl 
nil 



in a 



gh resolutio n HST image (|Chapman et al. ||2004 
Bower et al.||200l| . 



LiAB2 is associated with five if-band components. #1 
is a DRG and the brightest in if-band, but the photo- 
metric redshift is slightly lower than our current criteria. 
#2 is a DRG and the second luminous i^-selected object. 
There is a LBG at z = 3.091 in vicinity (indicated with 
a blue circle), 0".9 apart from #2. #3, #4 and #5 are 
newly detected objects. An X-ray source is detected near 
the position of the LBG, suggesting that there exists an 
AGN. MIPS 24 /jm source is not detected at the position 
of the if-band sources. 

LAB3 is associated with a single DRG, while two 
photo-z selected objects are located in vicinity in the 
east. This LAB is detected in X-ray, suggesting an AGN 
activity. 

LAB5 is also associated with a single DRG. This ob- 
ject is detected in X-ray, suggesting an AGN activity. 
The shape of Lya emission is elongated, or filamentary. 

LAB7 is associated with the three K-band compo- 
nents. Detect i on of a LBG at z — 3.098 was reported in 
Steidel et al. (2003), but they are separated in two ob- 
jects in our high resolution i^-band image (#1 and #3). 
#2 is also a LBG at z = 3.093. 

LAB8, which can be to be a part of LABI, has a LBG 
at z = 3.094. The if-band counterpart is newly detected. 

LAB11 includes three photo-z selected objects within 
the Lya halo. The spectroscopic redshift of the unre- 
solved optical counterpart is 3.0748. 

LAB12 has three photo-z selected objects. #1 is a 
DRG, and a LBG is marginally separated in two objects 
(#2 and #3). An X-ray source is detected at the position 
of #1. #1 is one of the most massive galaxies associated 
with LABs. 

LAB14 is associated with a singl e DRG, which is de - 



tected in X-ray and MIPS 24 ^m (Webb et al. 2009) 



One of the most massive galaxies associated with LABs. 

LAB16 has the two photo-z selected objects within 
the Lya halo. #2 is newly detected. 

LAB20 is associated with a single photo-z selected 
object, which is a LBG at z = 3.118. The deficit of Lya 
emission is seen at the position of the K-band source. 

LAB24 is associated with a iC-band source. Unfortu- 
nately, the color is not available as it is located at the 
edge of our images. 

LAB27 has the four photo-z selected sources. The 
shape of Lya emission is filamentary. 

LAB30 has the two photo-z selected sources. #1 is a 
LBG at z = 3.086. #2 is newly detected. 

LAB31 is associated with a single photo-z selected 
object, which is a LBG at z = 3.076. The shape of Lya 
emission seems round. 

Except for LABI, the most massive 4 LABs (LAB2, 
3, 12, 14) are detected in X-ray, which suggests that 
AGN activity. In addition, the luminous seven LABs 
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have DRGs in their vicinity, suggesting that the LABs 
are associated with dusty starburst galaxies. 

The detection of the multiple if -band components in 
LABs suggests that we are witnessing "on-going assem- 
bly" or " hierarchic al multiple merging" events of massive 
galaxy formation (White fc Rees 1978 Kauffmann fc 



White ||T993t|Cole et al.| I J0M[ |Meza let al.|M^|JNaab"eI 
al.|2007). The detailed discussion is given in section 5.2 



is 
ier 

12, 



4.2.2. LABs with no K -Counterparts 

Five LABs, LAB9, 19, 25, 26, and 35, which have 
no photo-z selected counterparts, are shown in Figfi"2") 
Those LABs have no luminous NIR and optical counter- 
parts in Lya halos, while LAB9, 19 and 35, have optical 
components within a halo which are classified as fore- 
grounds in our analysis. In addition, LAB9 and LAB26 
seem to be a close pair on the sky as shown in the right 
bottom panel in Fig[l2j while the object located between 
LAB9 and LAB26 is also supposed to be a foreground ob- 
ject in our analysis. Add itionally, the five L ABs are not 
dete cted either in X-ray ( |Geach et alT||2009[ ) or MIPS 24 
/zm ( Webb et al.| 2009). These LABs are supposed to 
be less massive objects with at least < 10 10 M . These 
objects ca n be predominantly e xcited by the cold gas 
accretion ( |Goerdt et "aL] ( |2010[ )). |Matsuda et al.| ( |2011 1 
suggested that the filamentary morphology of the nebu- 
lae may be the sign of cold accretion, although not all of 
the five LABs have such filamentary structures. 

Among them, only LAB35 has plausible counterparts 
in vicinity, which are one LBG and two DRGs. The 
redshift of the LBG is spectroscopically confirmed as 
z = 3.027, which is a bit outside of the protocluster sys- 
tem. In addition, one of the two DRGs, B in FigfLT 
classified as a foreground in our photo-z. On the ot 
hand, the other DRG, which is indicated as A in Fig 
has the photo-z in the range of z p hot = 2.6 — 3.6 an 
associated with an emission-line object. The question is 
whether the object A is related to LAB35 or not. The 
redshift of the object A (either the mission line compo- 
nent or the continuum) has not been determined spec- 
troscopically, while LAB35 is already known to be at 
z = 3.098. Their projected separation is 66 kpc if they 
are at the same redshift. The Lya emission of the two 
objects does not show any indication of their connec- 
tion; the emission at the object A and LAB35 are clearly 
snapped from each other and there is no counterpart be- 
tween them in any wavelength. 

5. DISCUSSION 

5.1. Massive galaxies in the SSA22 field 

In order to reveal whether massive galaxies have al- 
ready been formed in the high redshift protocluster, we 
investigate the distribution of the stellar mass of the 
if -band selected objects in the SSA22-M1 to M6 fields. 
The stellar mass is derived from the SED fitting of the 
BVRi'z'JHK photometric data. Spectral types, ages, 
and dust extinction are treated as free parameters. The 
redshifts are fi xed at the photometric redshifts described 
in section [3.1 1 while the spectroscopic redshifts are used 
for LBGs if available. The redshifts of the A"-selected 
objects in LABs described in section |4.2.1| are fixed at 
z = 3.1. W e used the template spect ra obtained from 
GALAXEV (iBruzual & Charlot||2003|) with star forma- 



tion histories with e- folding time r =0, 1, 9 Gyr. The 
best-fit SED is determined from the minimum \ 2 value. 
It should be noted that only the objects fitted with 
X 2 < 1 (85 % of the sample) are used for the following 
discussion. 

Fig[l3] shows the sky distribution of the if-selected 
galaxies with M* > 10 11 M . These red massive galaxies 
are supposed to be the representative population in the 
SSA22 protocluster. While the number density of the 
galaxies with M, > 10 11 M Q in GOODS-N is 0.16±0.04 
arcmin -2 , it was 0.37 arcmin -2 in SSA22, 2.3 times 
larger. 

In Fig |13| the surface number density of the massive 
galaxies is the highest near the peak of the LAE den- 
sity (cross) although a little bit shifted to the south. 
There are 14 galaxies within 4'x7' (« MOIRCS FoV) 
area, which gives the surface density of 0.5 arcmin -2 or 
3.1 times of the average of the GOODS-N field. 

On the other hand, we note that the number densi- 
ties of DRGs and massive galaxies in SSA22-M1 are the 
smallest among the observed six fields. This is interest- 
ing as the SSA22-M1 field is the densest region of LABs 
and the two largest LABs (LABI and LAB2) locate in 
the area. The field is ~10 Mpc (in comoving scale) away 
from the density peak of LAEs. 

Massive galaxies in protoclusters at z > 2 have been 
identifie d in some o ther fields mostly around the radio 
galaxies. Steidel et al.| ( J2005J ) identified a protocluster at 
z = 2.3 in the HS 1700+643 field, which contains a sig- 
nificant number of old and massive galaxies. From their 
spectroscopic sample in the ~ 70 arcmin 2 field, ther e are 
at least thr e e mem ber galaxies with > 10 11 M G 



©■ 



Ko- 



dama et al. (20071 observed the field around the high- 2 
radio galaxies at z = 2 — 3 with MOIRCS. They reported 
that there is an excess of galaxies with > 10 11 M on the 
bright end of the NIR color-magnitude sequence in the 
fields of MRC 1138-262 at z = 2.2 and USS 1558-003 
at z = 2.5 while there are few massive galaxies on the red 
sequence in the protoclusters at z = 2.9 and 3.1. For the 
MRC 1138—262 field, two of the massive red galaxies are 
identified as the protocluster m embers by MOIRCS spec- 



troscopy (jDohertyetal 
investigated the H5T IN 



2010b . Zirm et al. (2008) also 
CMOS images of the field and 



found the significant overdensity of the red faint galax- 
ies in the co lor-magnitude diagram. In addition, Hatch 



et al. (2009) derived the stellar mass of 19 galaxies in 
the Lya fialo of the radio galaxy MRC 1138—262, which 
is called as the Spiderweb galaxy (Miley et al. 2006), 
and found that there are no galaxies with > 10 11 M in 
the 0.3 arcmin 2 field except for the radio galaxy itself, 
which has the stellar mass of 10 12 M Q . Compared to 
these previous studies, thanks to our wide field coverage, 
our results clearly identified the overdensity of massive 
(> 10 11 M Q ) red galaxies in a protocluster at z = 3.1. 
Since this SSA22 protocluster is found as the largest high 
density region of star-forming galaxies known so far, the 
stella r m ass assembly likely to be more progressed. 

Fig[l4| shows the histogram of the stellar mass distri- 
bution for the photo-z selected galaxies. The number 
density in the SSA22 peak as well as that in GOODS- 
N is shown with the solid and dotted lines, respectively. 
Both the samples are selected with the same selection 
criterion with z p hot = 2.6 — 3.6 and K < 24.0. The turn 
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other redshift range (z p hot < 2.6 or Zphot > 3.6). The objects with numbers are adopted as if -band counterparts of the LABs. The blue 
circle in LAB2 is centered on the position of LBG M14 in Steidel et al. (2003). 



off of the distribution at the less mass end (log M*/M 
< 10.5) is due to the magnitude cutoff. The stellar mass 
distribution in GOODS-N without the cutoff is also in- 
dicated with the dotted thin line. It is found that the 
stellar mass distributions in both holds are not much dif- 
ferent, while the number density in SSA22 is larger than 
that in GOODS-N. The median values of the stellar mass 
with K < 24.0 in SSA22 and GOODS-N are 2.7 xlO 10 



M and 2.2 x 10 



10 



Mq, respectively. 



The excess in the surface number density of SSA22 to 
that of GOODS-N, which is regarded as the cluster com- 
ponent, is shown in the bottom panel of Fig 14j It might 
be interesting to compare this with the stellar mass dis- 
tributions in nearby typical rich clusters. As an example, 
the stellar mass distribution for the Coma cluster is also 
shown in the figure. We focus on the massive end (log 
M*/M Q > 10.5) of the mass function because of the mag- 
nitude cutoff mentioned above. The stellar mass function 
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Fig. 11. — continued. 



of the Coma cluster galaxies is derived from the .ff-band 



lumino sity function (LF) obtained by de Propris et al. 
(1998) (corresponds to M* K = -24.0, a = -0.78 for the 
bright end) and the typical mass to light ratio of local 
elliptical galaxies (M*/Lk — 0.6). The number density 
is averag ed in the observed area o f 29'. 2 x 22'. 5 (0.8 x 0.6 
~p98l 
b 



e 



the > 10 10,s M Q galaxies in SSA22 is equivalent to 
30-40 % of that in the Coma cluster. 

5.2. Stellar mass assembly in LABs 



Mpc) in de Propris et al 



physical density at 



and converted to the 
It should be noted that the 



weight for the brightest galaxy in the Coma cluster in 
fitting LF seems low (Fig. 2 in de Propris et al. 1998) 
and its contribution is not well reproduced by the ob- 
tained Schechter parameters. Therefore, we added it to 
the brightest bin manually. As shown in the figure, it 
is suggested that the equivalent of « 10 — 15% of the 
galaxies with > 10 105 M Q in the Coma cluster is al- 
ready seen in the SSA22 field. In practice, however, the 
central region of the Coma cluster is supposed to be ex- 
tended sever al times larger at z ~ 3. For example, |Suwa| 
et al. ( 2006 ) showed that a present-day cluster of galax- 
les with several Mpc can be extended to 20-40 Mpc in 
the comoving scale at z = 4 — 5 in their simulation. By 
assuming the spherical model for the nonlinear collapse 
of the overdense region, the surface number density at 
the epoch of turn-around is four times lower than that of 
the collapsed object. If the protocluster does not reach 
the maximum expansion, it may be slightly larger. As 
for crude estimation, the stellar mass distribution in the 
Coma clust er d ivided by three is shown in the bottom 
panel of Figjl4| The figure suggests that the number of 



As shown in section 4.2 most of the LABs are associ- 
ated with the -ftT-selected objects. The significant num- 
ber excess around LAB that these galaxies are formed 
in LABs. It is also interesting to see a notable frac- 
tion of them consists of multiple -ftf-selected components. 
The formation of early-type galaxies was generally un- 
derstood by either "Monolithic collapse" (Eggen,Lynden- 
Bell, Sandage 1962) or "(Major) Merger" (Searle & Zinn 
1978) scenario. In the realistic cosmological scenario, 
however, early-type galaxies may have experienced hier- 
archical gaseous multiple merging in their early evolu- 
tionary phase, which may be observed as the combina- 
tion of the extended gaseous halo and the multiple stellar 
components. 

Indeed, such phenomena is discussed using the cosmo- 
logical numeri cal simulations for the formation of early- 
type galaxies flMeza et al ||2003| |Naab et al.|2007| . The 
scenario can be called as "on-going assembly or "hierar- 
chical multiple merging" . The most outstanding exam- 
ples observed in our observational data are the two giant 
objects, LABI and LAB2, which have 6 and 5 X-selected 
objects, some of which are DRGs, within their Lya halo. 
As shown in the model simulations, these multiple com- 
ponents will finally merge into a single galaxy, which may 
evolve into a massive elliptical in the local universe. 

It is interesting to note that SSA22-M1 field where the 
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at z = 3.027 (see text). 
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Fig. 13. — The sky distribution of photo-z selected objects with 
Mi, > 10 11 Mq (filled circles, red). LABs are indicated with large 
open circles (blue). The lcr and la density levels of LAEs are 
shown with the black and cyan contours. The cross shows the 
density peak of LAEs. 



LABs are the most populate is a few Mpc (in physi- 
cal scale) separated from the center of the clustering of 
the if -selected galaxies or LAEs. The many (5/8) of 
the LABs in the SSA22-M1 field, LABI, LAB2, LAB7, 
LAB16, and LAB30, in fact have the multiple if -band 
components. While the galaxy formation in the central 
region of the cluster proceed earlier statistically, we may 
see more dynamically young phase of massive galaxy for- 
mation in the region a little separated from the center. 

We investigate the stellar mass of the photo-z selected 
objects associated with the LABs in order to investigate 
how much of stellar mass in these LABs have already 
been formed. The method is the same as that in the pre- 
vious section, while the redshift is fixed to be z = 3.1. 
The result is listed in Table 2J We here assume that all 
the broad-band emission is due to the stellar light and ig- 
nored the contribution of the nebula emission lines, which 
will be constrained by the future spectroscopic observa- 
tions. 

The stellar mass of the individual components of the 
LAB counterparts ranges from 3.8 x 10 9 M Q to 1.7 x 10 11 
Mq. Fig 15 also shows the distributions of the stellar 
mass of individual components and t he in tegrated stellar 
mass associated with the LABs. Figfl6| shows the color 
and the stellar mass of each individual component. LBGs 
are shown with the blue dots. As seen in the figure, low 
mass components in the LABs are dominated by LBGs. 
LAB7 and LAB11 have multiple LBGs with < 10 10 M , 
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text). The same but divided by a factor of 3 is indicated with the 
cyan dot-dashed line. 



which are close to each other. 

With respect to the integrated total stellar mass in 
each LAB, 70 % of the LABs (14/20) have > 10 10 M Q , 
and the median value is 3.1 x 10 10 M Q . The 20% of the 
LABs have the stellar mass M» > 10 11 M . Typical 
LBGs in the general field have M* ~ 2 X 10 10 M Q (Shap- 
ley et al. 2001), which is comparable or slightly smaller 
than the median stellar mass of the LABs of the current 
sample. 

The total stellar mass in each LAB versus the proper- 
ties of the Lya emission are plotted in Figjl7[ The hor- 
izontal axes show the integrated stellar mass of photo-z 
selected objects associated with each LAB, i.e., the ob- 
jects in Table[4] The vertical axes in Fig 17 (a), (b) and 
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(c) show the luminosity, size, and velocity width of the 
LABs, respectively. Note that there is a correlation be- 
tween the luminosity and the isopotal area of the Lya 
emission of the LABs (Matsuda et al. 2004) and these 
are not totally independent. In Fig|17|(a), or (b), we see 
that some massive galaxies with M* ~ 10 11 M Q are asso- 
ciated with the brighter and larger LABs (LABI, LAB2, 
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LAB3). As also suggested in Uchimoto et al. (2008), this 
implies that the origin of the Lya emission in these LABs 
is related to the previous (or recent) star-formation ac- 
tivity represented by the accumulated stellar mass. On 
the other hand, LAB12 and 14 have large stellar mass of 
M* ~ 10 11 M in spite of their low luminosity in Lya. 
Both LABs are associated with X-ray sources, and they 
are detected in 8 /zm and 24 ^m (Geach et al. 2009), sug- 
gesting that they suffer from dust obscuration. LAB 14 
is also detected as a SMG (Geach et al. 2005). In these 
cases, Lya emission is supposed to be suppressed by the 
dust. 



In Fig 17 (c), we can see that the total stellar mass and 



the velocity width of Lya emission for 9 LABs obtained 
by Matsuda et al. (2006, also private communication) 
have the notable correlation. Here we should note that 



the velocity widths are estimated by assuming a single 
Gaussian profile, and the stellar mass estimation includes 
the uncertainty of photo-z. Despite such uncertainty, the 
correlation suggests not only that the photo-z selected 
objects are indeed associated with the LABs but also 
that the estimated velocity width represents their phys- 
ical quantities. 

Although the c orrelation can be expla ined by some dif- 
ferent scenarios (Matsuda et al. 2006), we first discuss 
about the case in which we assume that the Lya gas 
clouds are gravitationally bounded. With the assump- 
tion, their dynamical mass can be evaluated from the 
velocity width and the radius which is represented by 
the extent of the Lya emission. Then, our result sug- 
gests that the dynamical mass of LABs is correlated with 
their stellar mass. It is known that the dynamical mass 
of local spheroidal galaxies at z < 1 correlates with the 
stellar mass by a f actor of ~ 5 — 10 (Drory et al. 2004 
Bundy et al.|2007 ). On the other hand, the stellar mass 
of most of our sample are in the range of 10 10 3 — 10 115 
Mm , and dyn amical mass of LABs estimated by Matsuda 
et al.| ( |2006[ ) ranges over lO 117 - 10 133 M . They are 
different by a factor of ~ 50 — 80. As the ratio seems 
significantly larger than the local value, the stellar mass 
assembly may be still in progress in LABs. 

Of course, LABs may not be virialized, and the virial 
radius may not be equal to the size of the Lya halo. 

As an alternative case, if we assume that the Lya emis- 
sion is due to the outflow from galactic superwinds, the 
characteristic timescales of star formation can be esti- 
mated by their spatia l extents and typical outflow speeds. 



Matsuda et al.| ( |2006[ ) evaluated the timescale of ~ 10 7 ~ 
yr tor the SSA22 LABs. Assuming the star formation 
rate of 10 3 M yr -1 as inferred from the stacked sub- 
mm flux, the estimated amount of the formed stellar 
mass, ~ 10 10 M©, is consistent with the timescale. In 
this case, it is suggested that the observed stellar compo- 
nents is formed from the star formation activities which 
cause the outflow. 

5.3. Dusty starbursts in the SSA22 field 

if-band selection is also useful to detect the highly 
dust obscured starbursting galaxies. They are important 
population which dominates the stellar mass of the pro- 
tocluster. We here focus on the dusty starbursts in our 
sample, namely, HEROs which are the object s wit h ex- 
tremely red color of J — K > 2.1 (see section [331). We 
also d i scuss for the AzTEC sources detected in I'lamural 
et al. ( |2009[ ), which are claimed to have the spatial cor- 
relation with Lya emitters. 

5.3.1. Overdensity of HEROs in the SSA22 field 

We detected 31 HEROs down to K = 24.0 in our ob- 
served field. The left and middle panels of Fig 18 show 
the sky distribution of the HEROs. The figure also shows 
the 2 a density contour of LAEs. It is found that the 
number density of HEROs in SSA22-M6 is the highest 
(12 HEROs) among the six fields in SSA22. In the right 
panel of FigfLS) the surface density of HEROs is shown 
as a function of a radial distance from LABs, LAEs, and 
AzTEC sources (see the next section) . The mean number 
density of 0.28 arcmin -2 in the SSA22 field is 2.3 times 
larger than that in GOODS-N (0.12 ± 0.03 arcmin^ 2 ) 



Assembly of Massive Galaxies in a High-z Protocluster 



15 




OLABs 
LAEs 



Q_ 



334.4 334.3 

R.A. (deg) 




o° 

334.5 334.4 334.3 

R.A. (deg) 



HEROs 



from LAB 1,2 

from LABs except LAB 1,2* 

from LAEs 

from AzTEC sources 

mean density 

GOODS-N 



>: 



±JL 



Radius (arcmin) 



Fig. 18. — (left) Sky distribution of HEROs. HEROs are indicated with filled circles (red). LABs and LAEs are indicated with larg e 
open circles (blue) and open squares (cyan), respectively. The contours show the high-density region of the LAEs ( |Hayashino et ah 2004). 
The inner contour shows the 2 a level. The cross is the position of the density peak of LAEs. (middle) AzTEC sources are indicated with 
solid open circles (black), (right) The surface density of HEROs in a circular ring of 0.1 arcmin width as a function of radial distance from 
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while those in SSA22-M6 is the 3.1 times larger. In addi- 
tion, the marginal excesses of HEROs are found around 
LABs/LAEs, while there are no HEROs around LABI 
and LAB 2. 

Out of the 23 HEROs detected in i?-band (see section 

33}, 11 have Zp hot = 2.6 - 3.6. The result of the SED 
fitting for the 11 HEROs shows that the stellar mass 
range from 2.7 x 10 10 M Q to 3.7 x 10 11 M , and the 
median is 1.6 xlO 11 M . The 55 % (6/11) of the HEROs 
are massive galaxies with > 10 11 M Q . It is likely that 
HEROs are among the most massive populations in the 
protocluster. 

Of all the 31 HEROs, two objects are associated with 
Lya emission, which are expected to be indeed located 
at z = 3.1. They are found in LAB 14 and the extended 
LAE in the vicinity of LAB35 (The object A; see section 
4.2.21. Both of them are located in the 2cr high density 



contour of LAEs. The object A in vicinity of LAB35 has 
J — K = 2.45 and z p hot = 3.06. The stellar mass of 
this HERO is 2.65 x 10 11 M Q , which is one of the most 
massive galaxies in this region. 

5.3.2. Correlation between K -selected objects and AzTEC 

sources 

The SSA22 protocluster is observ ed at 1100 iim usin g 
AzTEC camera mounted on ASTE ( |Tamura et al. 2009). 
They reported the concentration of su bmillimeter- bright 
galaxies (SMGs) in this region. Fig 19 shows the sky dis- 
tribution of the if-selected objects as well as the AzTEC 
sources in Tamura et al. (2009). The left panel indicates 
the distribution of the photo-z selected objects, and the 
middle one indicates that of DRGs. As shown in the 
right panel of FigjT9j the association between iif-selected 
objects an d the AzTEC so u rces is notable. 

Indeed, 



Tamura et al. 



(120101) reported the multi- 
wavelength identification of the brightest source, SSA22- 
AzTECl. The object is identified with the VLA radio 
source, whose counterpart is detected in K- and .ff-band 
images. On the other hand, no counterpart is detected 
in J-ban d and optical bands. The object is classified as 
a HERO. Tamura et al. (2010) estimated the redshift as 



•Zphot = 3.19 by their radio-submm photo-z. 



In the right panel of Fig[l8j we show the correlation 
between HEROs and the AzTEC sources. It is, how- 
ever, found to be dominated by the contribution of K- 
selected objects around SSA22-AzTEC20 as shown in 
FigJ20] They are supposed to be the counterparts of 
AzTEC20 and also in the "hierarchical multiple merg- 
ing" phase. It is interesting to note that significant Lya 
emission is not detected around the objects. 

Thus, the statistical evidence (clustering) as well as a 
few notable individual examples implies that at least a 
fraction of the AzTEC sources are indeed associated with 
the SSA22 protocluster. 

6. CONCLUSION 

We present the results of wide-field deep NIR imag- 
ing in the z = 3.1 protocluster region. The observed 
field is 112 arcsec 2 in the SSA22-Sbl field. We extracted 
the candidates of the protocluster members from the K- 
selected sample by using the multi-band photo-z selection 
as well as the color cut for DRGs, and examined their sky 
distributions. 

The surface number density of DRGs in our observed 
field shows an excess compared with those in the blank 
fields, which implies that progenitors of massive galaxies 
are growing in the protocluster at z — 3.1. Interest- 
ingly, the densest area corresponds to the density peak 
of LAEs, and the density is 2 times higher than that in 
the GOODS- North field, while most of the DRGs are not 
directly associated with the LAEs. On the other hand, 
the overdensity of DRGs in the SSA22-M1 field, where 
is characterized by the density excess of LABs, is not 
significant. 

We also investigated the near-infrared properties of 
LABs. if-band counterparts with z p hot — 3.1 are de- 
tected for 75 % (15/20) of the LABs within their Lya 
halo, and the 40 % (8/20) of LABs have multiple com- 
ponents. The stellar mass of 15 LABs range from 4 x 10 9 
M Q to 2 x 10 11 M Q . We also found that the stellar mass 
of LABs correlates with their luminosity, isophotal area, 
and the Lya velocity widths, implying that the physi- 
cal scale and the dynamical motion of Lya emission are 
closely related to their previous star-formation activities. 
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Fig. 20.— The if-selected objects around AzTEC20. The circles 
(black or red) indicate the photo-z selected objects with z p h t = 
2.6 — 3.6. DRGs are colored red. The two thick red circles are 
HEROs with Zphot = 2.6 — 3.6. The red square is a HERO with 
Zphot slightly lower than 2.6. The large circle shows the 2<r error 
range of the position of AzTEC20. The size of the image is 23". 4 
on a side. 



The results suggest that the formation of massive galax- 
ies is in progress in a large fraction of the LABs in the 
field. 

Overdensity of dusty starburst galaxies such as hyper 
extremely red objects (HEROs; J—Kab > 2.1) and plau- 
sible K-band counterparts of submillimeter sources are 
also found in this region. Highly dust-obscured galaxies 
which are expected to be growing up to massive galaxies 
may also be populated in such high density region of Lya 
galaxies. 

These results suggest that we are witnessing a proto- 
cluster in which the progenitors of massive galaxies are 
formed in the high-density region of star-forming galax- 
ies. 
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TABLE 1 
Summary of the observations 



Field 

(size) 



Center position 
R.A. Dec. 



Filter 



Date 
(HST) 



Seeing 
(arcsec) 


Exposure 
(hours) 


0.53 
0.48 
0.43 


0.9 
1.0 
1.5 


0.51 
0.47 
0.47 


0.9 
0.6 
1.1 


0.48 
0.59 
0.51 


2.0 
1.7 
1.2 


0.60 
0.50 
0.53 


2.0 
1.6 
1.5 


0.67 
0.46 
0.47 


2.1 
1.1 
1.3 


0.47 
0.54 
0.41 


1.2 

0.6 
1.0 



SSA22-M1 
(7' x 4') 



SSA22-M2 
(7' x 4') 

SSA22-M3 
(4' x 3'.5) 

SSA22-M4 
(4' x 3'.5) 

SSA22-M5 
(4' x 3'.5) 

SSA22-M6 
(7' x 4') 



22:17:33.9 00:12:14.0 



22:17:31.9 00:15:44.0 



22:17:16.6 +00:17:58.7 



22:17:52.8 +00:16:15.4 



22:17:07.4 +00:21:07.5 



22:17:28.7 +00:19:07.5 



J Jun. 15, 2005 

H Aug. 14, 2005, Jul. 22, 2006 

K s Jun. 14, 15, Aug. 14, 2005 
Jul. 22, 2006 

J Jul. 22, 2006 

H Jul. 22, 2006 

K s Jul. 22, 2006, Sep. 24, 2007 

J Sep. 22, 2007 

H Sep. 22, 23, 2007 

K s Sep. 22, 2007 

J Sep. 22, 23, 2007 

H Sep. 23, 2007 

K s Sep. 22, 23, 2007 

J Sep. 24, 2007 

H Sep. 24, 2007 

K s Sep. 23, 24, 2007 

J May 16, 17, 2008 

H May 17, 2008 

K s May 16, 2008 
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TABLE 2 
The limiting magnitude in each field 



Observed Area 


K a 


J b 


K b 


J c 


H c 


K c 


J d 


H d 


K d 




(l'V, 5a) 


(1.1" 


4>, 5cr) 


(1.4'V, 5 


*) 


(2.0'V. 5 


») 


SSA22-M1 


24.4 


24.6 


24.5 


24.4 


24.0 


24.1 


23.9 


23.5 


23.6 


SSA22-M2 


24.3 


24.6 


24.2 


24.2 


24.0 


23.9 


23.6 


23.5 


23.4 


SSA22-M3 


24.3 


25.0 


24.3 


24.7 


24.0 


24.1 


23.2 


23.6 


23.6 


SSA22-M4 


24.3 


- 


- 


24.3 


24.2 


23.8 


23.6 


23.5 


23.3 


SSA22-M5 


24.2 


- 


- 


24.2 


23.8 


23.9 


23.9 


23.3 


23.5 


SSA22-M6 


24.3 


24.7 


24.3 


24.3 


23.8 


24.0 


23.7 


23.4 


23.5 



a Measured in original images with seeing FWHM in Table IT] 
b Measured in smoothed images with 0".5 FWHM. 
c Measured in smoothed images with 0".7 FWHM. 
d Measured in smoothed images with 1".0 FWHM. 
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TABLE 3 
Near-infrared properties of LABs 



LAB 


K 


err 


J - K 


err 


^■phot 


^s pec 


Comment 


LABI 
















1 


22.02 


0.10 


1.24 


0.23 


3.38 






2 


22.97 


0.11 


1.47 


0.34 


3.08 


3.097 


DRG 


3 


23.34 


0.11 


1.15 


0.31 


3.18 


3.109 


LBG Cll a 


1 


23.69 


0.19 


>1.08 


- 


3.35 






5 


23.97 


0.22 


0.76 


0.51 


0.05 






6 


24.57 


0.23 


>0.78 


- 


2.96 






LAB2 
















1 


22.56 


0.10 


2.07 


0.38 


3.04 




DRG 


2 


22.91 


0.11 


1.68 


0.42 


3.00 


(3.091) 


DRG, vicinity of LBG M14 6 


3 


23.68 


0.20 


>1.03 


- 


3.25 






1 


24.30 


0.14 


0.41 


0.52 


2.98 






5 


24.51 


0.26 


0.75 


0.66 


3.65 






LAB3 
















1 


22.54 


0.12 


1.47 


0.29 


3.10 




DRG 


LAB5 
















1 


22.33 


0.10 


1.61 


0.19 


2.65 




DRG 


LAB7 
















1 


23.49 


0.12 


>1.26 


- 


3.19 


3.098 


DRG, LBG C6 


2 


23.65 


0.13 


0.19 


0.26 


2.65 


3.093 


LBG M4 


3 


23.97 


0.18 


0.71 


0.61 


3.35 


(3.098) 


(LBG C6), vicinity of 1 


LAB8 
















3 


24.35 


0.28 


0.66 


0.47 


3.35 


3.094 


LBG C15 


LAB 11 
















1 


23.79 


0.17 


0.65 


0.33 


2.79 


3.0748 


LBG C47 


2 


24.49 


0.23 


1.07 


0.47 


2.71 


(3.0748) 


LBG C47, vicinity of 1 


3 


24.94 


0.24 


0.59 


0.70 


3.35 


(3.0748) 


LBG C47, vicinity of 1 


LAB12 
















1 


22.11 


0.10 


2.14 


0.45 


2.12 




DRG 


2 


23.67 


0.21 


1.16 


0.35 


3.40 


3.094 


LBG M28 


3 


24.33 


0.18 > 


1.29 


- 


3.36 


(3.094) 


(LBG M28), vicinity of 2 


LAB14 
















1 


22.45 


0.11 


1.81 


0.34 


3.10 


3.089 


DRG 


LAB16 
















1 


22.99 


0.11 


0.92 


0.24 


2.71 






2 


24.22 


0.23 


>1.11 


- 


3.35 






LAB20 
















1 


23.30 


0.11 


0.93 


0.36 


2.58 


3.118 


LBG C12 


LAB24 
















1 


23.36 


0.13 


- 


- 


3.08 




DRG 


LAB27 
















1 


22.83 


0.12 


1.28 


0.26 


2.83 






2 


22.83 


0.12 


0.92 


0.34 


2.76 






3 


23.79 


0.13 


0.96 


0.40 


2.81 






4 


23.96 


0.17 


>0.89 


- 


3.42 






LAB30 
















1 


22.96 


0.11 


0.79 


0.20 


3.02 


3.086 


LBG D3 


2 


23.51 


0.19 


0.39 


0.31 


3.23 






LAB31 
















1 


23.20 


0.11 


1.36 


0.35 


2.87 


3.076 


LBG C4 


LAB35 
















A 


21.73 


0.10 


2.45 


0.26 


3.06 




DRG 


B 


23.04 


0.11 


1.85 


0.54 


0.84 




DRG 



ID of the LBG in Steidel et al. (2003) 
0".9 apart from the rest-frame UV source 
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TABLE 4 
The results of the SED fitting of 

PHOTO-Z SELECTED OBJECTS IN LABS 



LAB 


ID 


(1O 1O M ) 


x 2 


Mv 


LABI 


1 


14.49 


, +4.96 
-7.69 


0.70 


-23.01 




2 


9.15 


+5.02 
-5.83 


0.39 


-22.52 




3 


0.80 


+ 2.3!) 
-0.71 


0.28 


-21.95 




1 


2.41 


+3.41 
-2.15 


0.07 


-21.67 




5 


0.38 


+ 5.66 
-0.24 


0.10 


-21.68 




6 


1.14 


+ 1.69 
-1.12 


0.12 


-21.40 


LAB2 


1 


6.81 


+55.73 
—3.73 


0.03 


-22.91 




2 


10.92 


+ 10.25 
-10.05 


0.20 


-22.64 




3 


4.64 


+5.43 
-4.44 


0.05 


-21.80 




1 


1.44 


+ 1.52 
-1.33 


0.14 


-21.65 




5 


3.72 


+ 9.84 
-3.67 


0.01 


-21.31 


LAB3 


1 


9.26 


+7.30 

-4.28 


0.18 


-22.93 


LAB5 


1 


7.19 


+2.50 

— 2.5 5 


0.99 


-23.27 


LAB7 


1 


0.82 


+2.26 
-0.74 


0.12 


-21.72 




2 


0.75 


+2.47 
-0.68 


0.18 


-21.73 




3 


0.60 


+ 1.62 
-0.58 


0.06 


-21.31 


LAB8 


1 


0.42 


+ 2.32 
-0.41 


0.01 


-20.98 


LAB 11 


1 


0.42 


+ 1.35 
-0.39 


0.05 


-21.47 




2 


0.57 


+3.29 
—0.55 


0.07 


-21.30 




3 


0.58 


+ 1.65 
-0.52 


0.08 


-21.84 


LAB12 


1 


17.44 


+50.20 
-9.19 


1.15 


-23.02 




2 


2.07 


+ 1.93 
-1.86 


0.33 


-21.94 




3 


1.44 


+ 1.88 
-1.37 


0.31 


-21.65 


LAB14 


1 


16.95 


+5.19 
-9.91 


0.10 


-22.95 


LAB16 


1 


1.76 


+3.82 
-1.29 


0.55 


-22.48 




2 


0.95 


+ 1.50 
-0.92 


0.03 


-21.33 


LAB20 


1 


1.56 


+ 1.39 
-1.29 


0.17 


-21.86 


LAB24 


1 


2.85 


+2.78 
-1.69 


0.24 


-21.91 


LAB27 


1 


1.77 


+ 1.97 
-1.16 


0.85 


-22.56 




2 


3.32 


+2.08 
-2.85 


0.31 


-22.20 




3 


0.38 


+ 1.44 
-0.34 


0.25 


-21.56 




1 


1.77 


+3.01 
-1.70 


0.14 


-21.63 


LAB30 


1 


1.72 


+2.14 
— 1.02 


0.10 


-22.56 




2 


1.34 


+ 2.99 
-1.25 


0.19 


-21.80 


LAB31 


1 


2.56 


+ 1.97 
-1.59 


0.24 


-22.17 


(LAB35) 


A 


26.52 


+38.37 
-11.69 


0.12 


-23.74 



